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SUMMARY 

The  purpose  of  this  study  was  the  examination  of  an  empiri- 
cal method  of  obtaining  aircraft  transfer  functions  directly  from 
in-flight  records  without  recourse  to  lengthy  frequency  domain 
techniques. 

The  proposed  method  utilizes  a  precision  analog  computer 
and  plotter  to  generate  response  curves  matching  the  in-flight 
records* 

The  computer  set-up  used  for  simulation  is  so  arranged 
that  each  numerator  coefficient  of  the  transfer  function  is  directly 
related  to  a  potentiometer  setting.    At  the  most,  two  potentiometer 
settings  are  unknown,  since  the  first  and  last  coefficients  can  be 
obtained  from  the  initial  slope  and  steady  state  value  by  using 
Laplace  initial  and  final  value  theory.    In  the  case  of  longitudinal 
motion  the  period  of  the  phugoid  is  of  such  a  long  duration  as  to 
make  it  impractical  to  obtain  final  value  steady  state  conditions 
in  actual  flight  records.,    This  is  not  prohibitive  for  this  problem 
solution  method,  however,  since  one  coefficient  can  be  determined 
from  the  flight  record  initial  slope  or  suitable  accelerometer 
measurement  and  then  the  two  remaining  coefficients  are  de- 
termined by  adjusting  potentiometer  settings. 

Utilization  of  this  method  indicated  that  errors  would  be 
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of  the  order  of  five  per  cent  in  magnitude  with  the  combination 
of  computer  and  recorder  used.    It  is  felt  that  with  a  high  speed 
repetitive  computer  and  a  large  high-persistence  cathode-ray 
oscilloscope,  the  time  of  solution  and  degree  of  error  might  be 
reduced  if  sufficiently  accurate  flight  data  were  obtained 
initially. 

Further  study  is  indicated  for  inputs  differing  from  step 
functions  such  as  input  with  the  Laplace  transform: 

( '  -  s(sr-M) 

The  theory  for  this  case  is  mentioned  briefly  in  the  Proposed 
Technique  section* 
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INTRODUCTION 

Considerable  theoretical  work  has  been  done  with  frequency 
domain  techniques  as  applied  to  aircraft  dynamic  analysis,   whereas 
little  effort  is  indicated  in  the  field  of  transfer  function  determination 
directly  from  in-flight  transient  response  data.    The  purpose  of  this 
thesis  is  to  investigate  a  simplified  method  of  determining  transfer 
functions  by  the  latter  means. 

The  empirical  transfer  function  solutions  are  checked 
against  the  solutions  of  the  equations  of  motion  utilizing  stability 
derivatives  and  other  data  of  the  F-86D  aircraft.    This  aircraft 
was  chosen  since  it  is  of  relatively  recent  vintage  and  has  been 
in  existence  long  enough  to  supply  reliable  information  sufficient 
to  verify  the  techniques  developed  in  this  thesis,    Transfer 
function  poles  are  readily  determined  from  flight  test  records; 
however,  the  zeroes  require  tedious  computation  when  frequency 
domain  methods  are  used.    The  values  of  the  transfer  function 
poles  were  taken  from  analytical  computations,  however,  they 
were  checked  on  the  computer  for  validity  and  accuracy,    The 
previous  precaution  was  taken  in  order  to  narrow  the  field  of 
investigation  and  source  of  error  in  determining  the  zeroes  which 
were  the  crux  of  the  overall  problem „    It  was  subsequently  de- 
termined that  this  precaution  was  not  necessary  and  would  not 
be  required  in  an  actual  practical  situation. 


Two  basic  theorems  of  the  Laplace  transformation  (See 
Ref.   1)  constitute  an  important  part  of  the  development: 

(a)  Final  value  theorem  -  If  the  function  f(t)  and  its  first 
derivative  are  Laplace  transformable  and  if 

L  [f(t)]     r    F(s) 

then 

Limit  sF(s)    -    limit   f(t) 

s  -*•()  t— *»oo 

(b)  Initial  value  theorem  -  If  the  function  f(t)  and  its  first 
derivative  are  Laplace  transformable  and  if 

L[f(t)]    =    F(s) 

then 

Limit  sF(s)         limit  f(t) 

s  — >*co  t  — >-  0 

The  application  of  these  theorems  will  be  shown  in  con- 
junction with  measurements  taken  from  the  assumed  in-flight 
transient  response  curves. 

This  study  was  undertaken  byR.  Wo  Illgen  and  W.  P. 
Vosseler  under  the  guidance  and  direction  of  Professor  R.  M. 
Howe  as  a  joint  thesis  project  to  satisfy,  in  part,  the  require- 
ments for  a  Master  of  Science  Degree  in  Aeronautical  Engineering 
at  the  University  of  Michigan,   Ann  Arbor,  Michigan,  during  the 
period  September  1957  -  June  1958. 


SYMBOLS  AND  NOTATION 

CX  Anglo  of  attack,  degrees  or  radians 

b  Wing  span,  feet  or  inches 

c  Chord  length,  feet  or  inches 

CD  Drag  coefficient,   D/qS 

C,  Lift  coefficient,   L/qS 

C.^„~  Moment  coefficient  about  aerodynamic  center 

mac  J 

C^j  Normal  force  coefficient,  N/qS 

CN  Normal  force  coefficient  curve  slope  due  to 
angle  of  attack  dC^/dcx 

CM  Pitch  moment  coefficient,  m/qSd 

dCj, /dor  Slope  of  the  lift  curve 

&e  Elevator  deflection,  degrees  or  radians 

2 

g  Acceleration  of  gravity,  ft/sec 

h  Altitude,  feet 

L  Lift,  pounds 

m  Mass  of  airplane,   slugs 

M  Mach  number 

q  Dynamic  pressure,  >o  V  /2 

S  Wing  area,  square  feet 

V  Airspeed,  velocity,  knots,  mph,  ft/sec 

W  Weight  of  aircraft,  pounds 

/3  Angle  of  attack  in  yaw  plane,  degrees  or  radians 

CO  Natural  frequency,  rad/sec 


\  Damping  coefficient 

*y  Time  constant,   sec 


AIRCRAFT  AXIS  SYSTEM 


Axis 

X 

Force 

F 

X 

Moment 

L1 

Linear  Velocity 

u 

Linear  Acceleration 

u 

Angular  Displacement 

* 

Angular  Velocity 

P 

Angular  Acceleration 

• 

P 

Moment  of  Inertia 

I 
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yy 
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V     =    iU+  jV  +  kW         space  linear  velocity  vector 


SL     =    iPt-  jQ+  kR         space  angular  velocity  vector 


ILLUSTRATIONS 

Figure  Title 

1  Equipment 

2  Longitudinal  Computer  Diagram 

3  Longitudinal  Response  Curves 

4  Longitudinal  Simulator  Diagram 

5  Longitudinal  Simulator  Response  Curves 

6  Lateral  Computer  Diagram 

7  Lateral  Response  Curves 

8  Lateral  Simulator  Diagrams 

9  Lateral  Simulator  Response  Curves 


PROPOSED  TECHNIQUE 

THEORY 

Transfer  functions  can  be  obtained  from  the  equations  of 
motion  which  are  ratios  of  polynomials  of  the  general  form: 

M-   b^+b^B*-1* b0 

<f(s)        amsm  +   am_lSm^l  t       .......  ^ 

Multiplying  the  above  by  s  and  employing  the  Laplace 

initial  value  theorem  for  a  step  input: 

limit      s(sm-n)X   =    limit     dm""X(t) =■  -^—  <T 

s — >oo  t— *-o  dt  am 

The  slope  of  the  (m-n-1)  derivative  response  curve  is  then 

equal  to  bn^/am» 

The  final  value  theorem  for  a  step  input  gives: 

limit        sX    =    limit     X(t)    =  — °- — <f 
s  _^0  t — *.oo  a° 

The  first  and  last  coefficients  of  the  numerator  can  now 
be  easily  evaluated  since  the  denominator  coefficients  are  found 
by  the  simple  means  of  measuring  damping  and  frequency.    For 
n   =    2  there  are  (n-1)  unknown  coefficients  remaining. 

If  the  above  transfer  is  mechanized  (n-1)  potentiometers 
representing  the  unknown  coefficients  can  be  varied  so  as  to 
obtain  results  matching  the  in-flight  records.     This  is  not  overly 
difficult  to  accomplish  since  there  are  a  maximum  of  two 
potentiometer  settings  which  are  unknown  and  these  are  known 


at  least  to  an  order  of  magnitude.     Tins  would  apply  to  rigid- 
body  aircraft. 

A  suggested  method  of  procedure  for  a  step  function  in- 
put is  to  match  phase  fairly  closely  and  then  match  initial  slopes 
for  the  same  steady  state  value.    Iteration  brings  relatively 
rapid  convergence. 

For  the  case  where  the  input  transform  o  (s)  is         0 

s(sr   1) 

the  final  value  theorem  gives  results  identical  to  the  previous  case. 
The  initial  value  theorem,  however,  indicates  that 

limit  d—1  X(t) m       bn^T 

t    y  o  dt™-11-1  T  am 

where    /     is  the  time  constant  of  the  input. 

Only  examples  having  step  function  inputs  will  be  con- 
sidered in  this  thesis. 


EQUIPMENT 

The  computer  used  for  these  studies  was  a  20  amplifier 
electronic  differential  analyzer  designed  and  built  by  the  De- 
partment of  Aeronautical  Engineering,   University  of  Michigan. 
A  photograph  of  the  computer,  along  with  the  Variplotter,  is 
shown  in  Fig.   1*     This  computer  had  resistor  accuracy  of  .  02%, 

which  was  more  than  required  for  the  present  study.     The 

5 

computing  amplifiers  had  a  gain  of  approximately  10    and  were 

manually  balanced.,    Coefficient  potentiometers  were  set  to 
an  accuracy  of  0.  1%. 

The  Pace  Associates  Model  HOOD  Variplotter  used  to 
plot  the  analogue  computer  outputs  is  a  table- top  size,  self- 
balancing,  potentiometer  type  of  recorder.    It  records  two 
variables  simultaneously  in  the  form  of  a  rectangular  coordinate 
graph.     The  static  accuracy  of  both  the  pen  and  arm  is  .  075  per 
cent.     The  dynamic  accuracy  for  both  is  .  2%  for  plotting  speeds 
up  to  15  inches  per  second. 
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EXAMPLES 

Examples  of  the  determination  of  zeroes  for  both  longi- 
tudinal and  lateral  functions  will  be  given. 
I.     Longitudinal  Transfer  Function  Zeroes 

Since  no  suitable  flight  records  were  available,  flight 
data  was  assumed  from  solving  the  aircraft  equations  on  an 
electronic  differential  analyzer. 

Under  conditions  of  perturbations  from  steady  level 
flight  at  Mach  .  8  and  30,  000  feet  the  following  equations 
were  utilized  based  on  the  airplane  data  in  Appendix  A. 

ii    =    -.  00664u  +   .  0497w  -  31.  9      -  32.  2  e 

q 

w=  -,0407u  -  l,018wi-790.3q  -  1.  29  9  -  1.  55  ff 
q  =  .000475U  -  . 0119w  -  , 737Q         -  . 413  S 

e  -  q 

Scaling  the  set  of  equations  for  use  on  the  analog  computer 
gives 

u  -    -.  00664u  +  .  00497(10w)  -  .  319(100q)  -  .  322(100  6) 

w  =    -.0407u    -    .1018(10w)+ 7.  903(100q)  -  .0129(1006)  -.155(10$) 

100  q     =    .  0476u  -  .  119(10w)  -„  737(l00ci)  -.  0413(10$) 
100  6     z    lOOq 

These  equations  were  then  mechanized  as  shown  in  Fig„  2 
(see  Appendix  B  for  all  remaining  figures).     The  following  scaling 
was  used: 

1  fps    =    .  1  volts 
1  rad   =    .  1  volts 
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1  deg  =    .  1  volts 

1  rad/sec     =    .  1  volts 


The  exact  transfer  function  of  Yqs-    was  found  from 


this  set  of  equations  to  be 
Y    =    °/s 


23.  6(s2   +    •  983  s  +  .  00855) 
(s     +-   1,762s3  +     10,16s2    +   o  0195s  +.0311) 


Since  there  is  very  little  coupling  between  the  short  period  and 
the  phugoid  motion  it  was  considered  possible  that  the  poles  and 
zeroes  of  the  above  transfer  would  be  the  same  as  the  poles  and 
zeroes  of  the  uncoupled  transfer  functions. 

The  approximate  formula  for  the  short  period  transfer 
function  as  based  on  ReL  2  is 


Kg  ^(1.032  8  +    1) 
Vfe    z         (.0978  s*    +  .1885  s    +  1) 


The  approximate  formula  for  the  phugoid  transfer  function 
as  obtained  from  the  same  source  is 
KQSe  (127  s   +  1) 

6/cT,  = g 

e  (305.5  s     +  2.405  s    +-  1) 

The  product  of  the  above  two  expressions  does  indeed  give 

a  very  close  approximation  to  Ye  r    which  makes  a  lengthy 

°  e 

Newton  Raphson  approach  (see  Ref.  3)  or  other  factorization 
technique  unnecessary. 

The  parameters  of  the  approximate  transfer  function 
are  then  as  follows: 
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/  »        =    3.  22  radians 
w  ns 

<°np 

=    .0572  radians 

Ts      =    1.862  seconds 

TP 

=    109. 5  seconds 

!fs    =  °304 

rP 

=    .0685 

These  values  were  used  tc 

>  assure  that  the  response 

curves  obtained  from  the  mechanization  of  these  curves  were 

indeed  correct.     The  computer  response  curves  of  Fig.  3  had 

the  following  parameters 

CO        -    3. 22  radians  60n„    =    ,0572  radians 

ns  np 

Tc      =2  seconds  T„       =    112  seconds 

to  p 


P    •    •« 


h    -  -3  f, 

The  above  parameters  give  the  poles  immediately  in  the 


canonical  form0 

The  determination  of  zeroes  for  the  transfer  function 
Y  ©  r    was  considered  next. 

In  the  interest  of  a  simpler  circuit,  the  constant  term 
(.  0311)  in  the  denominator  was  assumed  to  be  negligible    An 
1  s*  was  then  factored  out  of  the  remaining  terms  giving  the  transfer 
function  of  Y 


q<f  e 


-23. 6(s2    -t  -983s    Hr  ,00855) 


^e  (s6  -»■    1.762s2   +  10. 16s  ■*-  .0195) 

which  has  the  dimensions  of  degree  per  degree  or  radian  per 
radian. 

The  denominator  is  now  a  cubic  with  coefficients  which 
can  be  equated  to  potentiometer  settings  and  amplifier  gains    as 
follows  for  the  longitudinal  simulator  computer  circuit  of  Fig.  4: 
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W1K1K2    z    1>762'  Kl    -    lt     K2    =    10,      1    =    -1762 

CX2K1K2K3    -    10.16,  K3    =    2,    0<2    =    .508 

OCjKjKgKgl&j    =    .0915,        (ocK)7     r    .  1,  0(3    r    .010 

Inasmuch  as  the  steady  state  of  the  phugoid  is  difficult 

to  obtain  in  flight  tests,  it  was  decided  that  a  more  practical 

method  of  procedure  would  be  to  determine  the  potentiometer 

setting  affecting  the  initial  slope  by  utilizing  flight  records  and 

to  vary  the  remaining  two  potentiometers. 

The  initial  slope  potentiometer  setting  corresponding  to 

2 
the  coefficient  of  the  s    term  in  the  numerator  was  found  in  the 

following  manner. 

The  initial  slope  is  found  from  Fig.  3a  or  3b.    Recalling 

that 

i</<T     =    b2/a3 
we  have  from  Fig.  3 

b2   =    .  4  rad/sec  /degree. 

This  was  used  to  obtain  the  potentiometer  setting  as 
shown  in  Appendix  C, 

As  has  been  indicated  before,   steady  state  for  the  phugoid 
takes  too  long  to  obtain  and  a  cut  and  try  approach  is  more 
expedient.    This  leaves  potentiometers  five  and  six  which  must 
be  adjusted  in  the  computer  setup  of  Fig.  4.    Various  combinations 
over  the  full  range  of  each  potentiometer  were  tried.    At  each 
setting  an  output  trace  of  q  for  a  1°  elevator  step  was  compared 
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with  the  assumed  flight  record.    Noting  the  effects  and  trends  of 
these  combined  settings  it  was  possible  to  come  within  5%  of 
what  the  actual  required  settings  should  be  after  only  about  18 
trials.    Increased  magnification  using  specialized  high-persistence 
cathode  ray  oscilloscopes  should  increase  accuracy  of  such 
measurements. 

Fig.  5  illustrates  the  technique  in  a  general  way  uti- 
lizing the  following  sequence  of  settings: 

Trial  Potentiometer 

5  6 

A  .5.5 

B  o75  .5 

C  1.0  .5 

D  1.0  .0 

It  is  to  be  noted  that  the  initial  slope  remains  unaffected 
as  might  be  expected.     The  latter  condition  suggests  that  the 
potentiometer  4  setting  could  also  have  been  obtained  by  trial 
and  error  matching  instead  of  from  initial  slope  measurement. 
Variation  of  potentiometer  5  affects  the  time  and  magnitude  of 
the  first  peak  as  well  as  setting  the  continued  level  of  the  entire 
curve.    Variation  of  potentiometer  6  has  a  negligible  effect  on 
the  initial  peak,  however,  it  raises  or  lowers  the  tail  of  the 
curve  thus  influencing  the  final  or  steady  state  value- 
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II.     Lateral  Transfer  Function  Zeroes 

Again,  since  no  suitable  flight  records  were  available, 
the  flight  response  curves  were  obtained  by  solving  the  air- 
craft equations  on  an  electronic  differential  analyzer. 

In  a  manner  similar  to  the  longitudinal  case,  the  air- 
plane data  in  Appendix  A  was  used  to  obtain 
v  =  -.  179v  +  31.  8p  -  796r  +  32.  2  <j> 

p=  -.435v  -  3.6p+  1.295r  +  ,313r  -  1.015cT    +  .164  J1 

a  r 

r  =  .  0121v  -  .  0389p  -  .  388r  +  .  0765p  -  .  0146  <T  -  .  0838  cT 

a        r 

fi  =  p  +  .  04r 

Scaling  the  set  of  equations  for  use  on  the  analog 

computer  gives: 

v  =  -.  179v   +■  .  318(100p)  -  7.  96(100r)  -t  .  322(100^) 

lOOp  =  -.  0435v-  3.  6(100p)  +   1.295(100r)  +   .  313(100r)  -  101.  5^ 

t  16.4J- 
r 

lOOr  =  1.  21v  -  .  0389(100p)  -  .  388(100r)    t  .  0765(100p)  -  1.  46  <T 


-8.38<f 


r 


100#  =  lOOp  +  .04(100r) 


These  equations  were  then  mechanized  as  shown  in  Fig.  6. 
The  scaling  was  as  follows: 
1  fps  =  o  1  volt 
1  rad  =  o  1  volt 
1  deg  =  .  1  volt 
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A  static  check  was  performed  by  observing  the  voltage 
on  inverters  1,  6,   and  11  for  various  initial  condition  voltages. 
The  response  curves  are  shown  in  Fig.  7. 

The  important  transfer  functions  for  this  set  of  equations 

were  found  to  be  as  follows: 

Y<A£      -    -.138(s2  -  2.45s  -  9.86) 

D(s) 

Y6  J"      -    1  051(s2  -h   .582s  +  10.1) 
a  D(s) 

Y  r         =  -62.7(s2  +    4.06s  -  .0296) 

V  r  D(8) 

Y  r       =    -42(s2  -  .872s  -,325) 

V  a  D(s) 

Y  r      =    .Q73(s3-t-  4.49s2-h  1.45s  -r-,752) 


rdr 


D(s) 


Y    r      =    .0945(s3   -I-  ,321s2  +  4.48s  ±  4.51) 

D(s) 

where 

D(s)    =    s4  ±   4ol6s3  +  10.  6s2+    38.5s  -  .01975. 

Since  the  transfer  functions  for  uncoupled  rolling  motion, 
spiral  motion,   and  the  dutch  roll  motion  promised  to  be  useful  in 
finding  the  poles  and  zeroes  of  the  above  transfer  functions  these 
modes  of  motion  were  investigated.    For  uncoupled  rolling  motion 
we  assumed  that  r  r/3  =  (f    -    °-    This  gives 

p   =  -3.6p  -  1.015cfa 
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or 

ft**-- >** • 

s(.  287s    +    1) 

For  the  dutch  roll  motion  we  assumed  that   p   =  (p  -  cT    =    °« 

This  gives 

v    =    -.  179v  -  796r 

r   =    .0121v  -  ,388r  -  . 0838  cT 

wr 


or 


.  0838  (s      t     .  179) 


/Or    r         Is2    -t-     .567s    t9.72) 

For  spiral  motion  we  assumed  that  v    =    p=r   =   /3=J^ 

<-T     =    0.    This  gives 
^  r 

0  =    -,179v  +  31.  8p-  796r  +  32.  2   <f> 

0  =    -4o35v  -  360p   t  129.  5r 

0  =    1.21v  -  3.89p  -  38.  8r 

0  =  <b    -p  -  .  04r 
If  we  assume  in  addition  that  v  =    p  =    0,  we  get  a  spiral  mode 
time  constant  of  1875  seconds. 

With  these  as  approximations,  synthetic  division  coupled 
with  the  Newton-Raphson  iteration  technique  for  determining  roots 
was  used  to  give 

Y(t>Sr    ~    ~'138<s  "  3.267)  (s  +  3.022) 

D(s) 

Y<*  f  a    =    1.051(s2  »   .582s  +  10.1) 
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Y  r      =    -62.7(s  J-  4.06)  3 
V     r  D(s) 

Y  f      =    -41>(s-  1.15)  (3    f  .282) 

D(s) 

Y  J"      r      Q73(s  +   4  1)  (s2  +  ,32s  j-  ,2) 
r     r  D(s) 

Y  <f     =      0945(s   t  .91)  (s2  -  .5793  f  5) 

a  D(s) 

where 

D(s)    =    (s  -  .000513)  (s  -r   3.95)  (s   +  .2255s  +    9.742) 


The  parameters  associated  with  the  denominator  of  the 
various  uncoupled  transfer  functions  are  then  as  follows: 
Tr   =    .  287  seconds  -  the  time  constant  of  the  uncoupled  rolling 

motion* 
T     =    1875  seconds  -  the  time  constant  of  the  spiral  mode 

ndr    -    • 324  radians  Per  second  -  dutch  roll 


J 


dr 


=    .875 


The  parameters  associated  with  the  denominator  of  the 

transfer  function  are  found  to  be: 

T     =    .  253  seconds 
r 

TQ    =    1950  seconds 
^ndr    =    .324 


875 


J^dr 

These  parameters  give  the  poles  immediately  in  the 
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canonical  form. 

The  parameters  of  the  simulator  circuit  shown  in  Fig.  8 
are  determined  by  setting  the  coefficients  equal  to  the  potentiometer 
settings  and  amplifier  gains  as  follows: 

oCjKjKg    =    4.16  Kj    =    1,  K2   =    10  «1    =    .416 

e^K^    -    10.6  K3    *    2,  <*2    -    .53 

«3K1K2K3K?    r    38.4      K?=    2,  <X  3    =    .  96 

since 

D<s)    =    s(s3  +   4.16s2  +  10.6s  +38.4) 

Two  transfer  functions,  YqJj-  and  Y •  rr  were  investigated, 
as  being  representative.    The  value  of  the  potentiometer  setting 
corresponding  to  the  highest  degree  's*  term  in  the  numerator  was 
found  as  in  the  longitudinal  case.    See  Appendix  C  for  details.    The 
remaining  potentiometers  were  adjusted  giving  the  plots  and  variations 
shown  in  Fig.  9, 
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RESULTS  AND  CONCLUSIONS 

By  determining  potentiometer  settings  corresponding  to 
initial  slope  and  also  the  steady  state  value,   where  available, 
it  was  possible  to  obtain  accuracy  in  matching  the  curves  such 
that  the  zeroes  were  determined  with  errors  of  less  than  5%, 

In  the  case  of  q/£e  ,  where  two  potentiometer  settings 
were  unknown,  only  18  trials  were  necessary  to  obtain  this 
accuracy. 

See  Fig.  5  and  Fig.  9  to  observe  the  variation  in  re- 
sponse curves  for  variation  in  numerator  potentiometer  settings. 
It  is  concluded  that  the  best  general  procedure  is  as  follows: 

1)  Obtain  estimates  of  poles  and  zeroes  from  pre-flight 
data. 

2)  Obtain  poles  from  flight  records  involving  step  inputs 
and  coefficients  of  first  and  last  terms  in  numerator 
corresponding  to  initial  slope  and  steady  state,  if 
possible.    Then  compare  with  previous  estimates. 

Limit      dm-nX/dtm-n=bn<r/am 

tr-*.    0 

Limit      X(t)  =  bQS/a0 
t — ►-co 

3)  Adjust  potentiometers  to  correspond  to  the  known 
parameters. 

4)  Set  estimates  in  those  potentiometers  which  must 
be  adjusted. 


21 


5)  Make  preliminary  adjustments  to  observe  which 
potentiometer  has  more  effect  in  the  case  where 
two  potentiometer  settings  are  unknown, 

6)  Adjust  the  more  critical  to  give  loose  matching, 
observing  trend  of  frequency,  phase  shift,  peak, 
and  steady  state  values. 

7)  Adjust  the  remaining  potentiometer  for  a  closer 
match  while  observing  trends. 

8)  Return  to  the  first  potentiometer  and  repeat  ad- 
justment. 

9)  Continue  iterating  for  as  close  a  match  as  desired. 
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APPENDIX  A 
AIRCRAFT  DATA  AND  STABILITY  DERIVATIVES 

The  following  data  is  given  in  references  4  and  5  for 
the  F-86D  aircraft  at  30,  000  feet,  Mach  0.8,  weight  15,  000  lbs. 


S  =  288  ft2 

1^  =  8625  slug  ft2 

Izz  =  35,  300  slug  ft2 

V    =  796  ft/sec 
P 

C*=  2.3°  -  0.04  rad 
/O  =  .  00089  slug/ft3 


CL  = 

0.185 

c 

MQ 

-  -6.5 

CL« 

=  4.  62/rj 

Cn 

=  .031 

l 

CM« 

=  -0.424 

CM± 

=  -1.5 

CD« 

=  .143 

cD- 

.018 

CM  I*    ~  ~-  0185/degree 
6  e 

CT    ~     =0. 0089/degree 
"  o  e 


c  =  8.09  ft 


Iyy  =  29400  slug  ft2 
Izx=  2,700  slug  ft2 
CT      =  -.82/radian 


b=  37,12  ft 


%  = 

C*T               - 

°YR- 

Y<Ta 

C„,    - 

V 

V 

-0. 10/radian 

V 

.117 

LP 

-.433 

C        = 
LR 

.15 

CNP  = 

-.01 

°Nr  = 

-.2 

°  a 

=  -. 0029/degree 

cL<rr' 

=  . 00043/ degree 

CN  c    =  "•  000055/degree 

CN<fr=-001 
CXT       =  .  117 

CT       =  -.10 

L/3 
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Fig.  2 
Cumputer  Circuit  for  Solving  the  Longitudinal  Flight  Equations 
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APPENDIX  C 


Determination  of  Potentiometer  Settings 

For  the  longitudinal  case  we  found  that 
b2  =  .  4  rad/secvdegree. 
from  Fig.  3  a  or  3b. 

The  scale  that  we  wish  the  plotter  to  have  is 
1  inch  -  .  1  rad/sec 
with  the  pen  scale  factor  set  to  10  volts  per  inch. 

For  a  one  degree  step  input  the  initial  q  should  be  equal 
to  40  volts,  since 

40  volts  =  .  4  rad/sec  . 


The  equation  relating  the  40  volts  to  the  parameters  of 
sui 

St 


the  circuit  can  be  determined  by  reference  to  the  figure  below„ 

Z_ 

DCS) 


IJ 


■[2/D(s)](^4K;(fe(sVK3K7^Sl 


Limit  sX(s)  =  qQ  =  (2/ag)  (o<4Kfc  J^/KgKy)  =  40  volts 
s >  00 

For  S     a  2  volts 

*4K6  =  2 

In  brief  then  for  the  initial  slope 

(K/a3)  (Of4K^  <Te/K3K7)  =  (K^/iy)  (bg/ag) 

and  for  steady  state 

(K/aQ)  («6k^  <Te)  =  (Kvi/Kri)  (b</ao) 


C-2 


where 

K  =  denominator  multiplying  factor 

K   .  -  the  pen  scale  factor  on  the  plotter  (in  the  above 

example  it  was  ten  volts  per  inch. ) 
Kri  z  the  number  of  radians  /second  or  the  number  of 

radians/second  /second  (in  the  above  example 

it  was  .  1) 


thesl3 

Determination  of  aircraft  transfer  func 
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